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Abstract 

Introduction: Visual but subjective reading of longitudinal relaxation time (T,) weighted magnetic resonance images is 
commonly used for the detection of brain pathologies. For this non-quantitative measure, diagnostic quality depends on 
hardware configuration, imaging parameters, radio frequency transmission field (S, + ) uniformity, as well as observer 
experience. Parametric quantification of the tissue relaxation parameter offsets the propensity for these effects, but is 
typically time consuming. For this reason, this study examines the feasibility of rapid 2D quantification using a variable 
flip angles (VFA) approach at magnetic field strengths of 1.5 Tesla, 3 Tesla, and 7 Tesla. These efforts include validation in 
phantom experiments and application for brain T, mapping. 

Methods: T q quantification included simulations of the Bloch equations to correct for slice profile imperfections, and a 
correction for 8, + . Fast gradient echo acquisitions were conducted using three adjusted flip angles for the proposed T, 
quantification approach that was benchmarked against slice profile uncorrected 2D VFA and an inversion-recovery spin- 
echo based reference method. Brain T, mapping was performed in six healthy subjects, one multiple sclerosis patient, and 
one stroke patient. 

Results: Phantom experiments showed a mean T n estimation error of (-63±1.5)% for slice profile uncorrected 2D VFA and 
(0.2±1.4)% for the proposed approach compared to the reference method. Scan time for single slice T, mapping including 
6, + mapping could be reduced to 5 seconds using an in-plane resolution of (2 x2) mm 2 , which equals a scan time reduction 
of more than 99% compared to the reference method. 

Conclusion: Our results demonstrate that rapid 2D T, quantification using a variable flip angle approach is feasible at 1.5T/ 
3T/7T. It represents a valuable alternative for rapid T q mapping due to the gain in speed versus conventional approaches. 
This progress may serve to enhance the capabilities of parametric MR based lesion detection and brain tissue 
characterization. 
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Introduction 

Magnetic resonance imaging (MRI) offers capabilities for non- 
invasive tissue characterization for a broad range of MRI 
applications including neuroinflammatory diseases and stroke 
[1,2,3,4,5]. At high and ultrahigh magnetic field strengths, brain 



parenchyma can be depicted with higher spatial resolution, which 
improves morphological conspicuity and facilitates improved 
characterization of multiple sclerosis (MS) plaques [6,7,8], 
enhances differential diagnosis of orphan neuorinflammatory 
diseases [9,10,11] and stroke [12]. Visual assessment of T! 
weighted techniques [13,14,15] is a common standard in today's 
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clinical brain imaging practice, although the diagnostic efficacy of 
cerebral lesion detection depends on the severity of focal or 
regional (patho-)physiological changes, on the imaging technique 
and pulse sequence design/ timing used, on the impact of 
hardware configuration, as well as on the observer experience. 
Notwithstanding the ubiquity and success of clinical T] weighted 
imaging, quantification of brain parenchymal Tj values is of great 
research interest and of substantial clinical relevance [16,17,18] 
but comes with the caveat of typically being rather time 
consuming. 

T[ quantification techniques using variable flip angles (VFA), 
such as driven equilibrium single pulse observation of Tj (DESPOT 1, 
[19]) have been proposed for fast three dimensional Tj mapping of 
the brain. DESPOT 1 includes at least two radio frequency (RF) 
spoiled gradient echo (FLASH [20]) measurements using different 
flip angles from which T! maps with full brain coverage can be 
achieved in approximately 10 minutes scan time [21]. Not all 
clinical indications necessarily require whole brain coverage, but 
might rather benefit from fast 2D zonal imaging covering 
particular brain regions or lesions using a high in-plane spatial 
resolution together with a limited number of slices. Consequently, 
targeting rapid 2D mapping approaches is conceptually appealing 
for clinical T! quantification of brain tissue. However, commonly 
used non-ideal radiofrequency (RF) pulses evoke an inhomoge- 
neous slice excitation leading to deformed slice profiles [22]. Slice 
profile deformation presents an extra challenge for VFA based T! 
mapping and bears the potential to deem T[ quantification 
inaccurate. Recognizing the constraints of conventional 2D VFA 
and the opportunities of 2D Ti mapping, this study examines the 
feasibility and fidelity of rapid slice profile corrected 2D VFA Ti 
mapping. To meet this goal, numerical simulations of the Bloch 
equations are used, which account for non-ideal RF pulse shapes 
and RF transmission field (Bj^) non-uniformities. The applicability 
of this approach is demonstrated in phantom experiments at 1.5T, 
3T and 7T. In vivo feasibility studies including healthy subjects, MS 
and stroke patients are conducted as a precursor to a broader 
clinical study. The merits and limitations of the proposed 2D VFA 
T! mapping variant are discussed and implications for clinical 
imaging are considered. 

Methods 

Theory 

Fast radio frequency (RF) spoiled gradient echo (FLASH) 
measurements using at least two flip angles allow for the 
quantification of the longitudinal relaxation time T! [19]. The 
commonly accepted equation that governs FLASH signal intensity 
S(oc) in the steady state is [23] 

5(a) = M 0 (1 - E 0 Si " ( °? e ~ffi, withfr = e "ff , (Eq.l) 
1— £iCOs(a) 

where M 0 is the proton density, TR is the repetition time, TE is 
the echo time, a is the flip angle, and T 2 * is the effective 
transversal relaxation time. Rearrangement of equation 1 yields 
[24] 

=£ > -EOe'^ (Eq. 2) 

sin (a) tan (a) 

Arranging S(oc)/sin(oc) over S(oc)/tan(a) allows for extraction of 
Ti from a linear fit T! = -TR/ln(m), where m is the slope between 
measurement points. 



Bloch equation simulations 

Equation 1 and 2 hold true for excitations that exhibit a uniform 
flip angle over the entire target area. In 2D acquisitions, however, 
non-ideal RF pulse shapes lead to deformations of the slice profile 
altering the resulting signal intensities. To assess the extent of slice 
profile deformation and its impact on variable flip angle Ti 
quantification, a Matlab (TheMathworks, Natick, MA, USA) 
environment was developed allowing for simulations according to 
the Bloch equations and emulation of a FLASH sequence. A 
bandwidth truncated sine RF pulse (bandwidth = 1 kHz, dura- 
tion = 2 ms, time-bandwidth-product = 2) used for FLASH imag- 
ing was extracted from the sequence development environment of 
the "integrated development environment for (MR) applications" 
(IDEA, Siemens Healthcare, Erlangen, Germany) and fed into the 
simulation. One-dimensional magnetization vector profiles along 
the slice selection direction were generated in the steady state for 
given flip angles. 

Ti quantification: Simulated versus experimental signals 

Using mathematical integration of the simulated slice profiles 
for different flip angles, theoretically achievable signal intensities 
were calculated. An iterative least squares minimization function 
matched these theoretical signal intensities to measured FLASH 
MR signal intensities by varying Tj between 1 ms and 10000 ms 
starting at T[ = 1000 ms and adapting the scaling factor M () . Pixel- 
by-pixel application of this method facilitated Tj mapping. The 
simulation environment allowed inclusion of afore acquired Bj + 
data that were derived from RF transmission field mapping. 

Three flip angles approach 

Further to the conventional approach that uses two flip angles 
we propose the application of three flip angles for VFA based Tj 
quantification. The first flip angle C*i was chosen to achieve the 
maximum signal intensity determined by Bloch simulations using 
the specific target T i . To determine the remaining two flip angles, 
oc 2 and oc 3 were varied from 1° to 90° including all possible 
combinations. Gaussian noise was added to the corresponding 
signal and T] evaluations using the proposed approach were 
performed for all flip angle combinations. Flip angle sets resulting 
in minimal Tj differences between the theoretical T, value and the 
estimated T 1 value were regarded as adjusted flip angles and were 
henceforth used as nominal flip angles for the phantom 
experiments and for the in vivo studies. 

In the first step, these one-time calculations were performed for 
Tj of a manganese chloride doped water phantom (Ti = 1000 ms) 
described below. In the second step, approximated mean T] values 
found in the literature for gray and white matter at different field 
strengths (T 1 = 950 ms at 1.5T [25,26,27,28], T! = 1250 ms at 3T 
[25,27,29,30,31], and Ti = 1650 ms at 7T [27]) were targeted. 

Monte Carlo simulations 

To assess T[ dependent accuracy (defined as deviation of mean 
Tj value from nominal Tj value) and precision index (defined as 
standard deviation of Tj divided by mean T t ) of the proposed 
mapping method and to compare the three flip angles approach 
with the two flip angles approach, Monte Carlo simulations of the 
proposed corrected 2D VFA fitting algorithm were performed for 
T] values of 1 ms to 3000 ms using a 1 ms increment. Gaussian 
noise was added to the signal intensities so that a SNR margin of 
100 between the highest of all three signals and the standard 
deviation of the noise was achieved. 10000 simulation experiments 
were conducted for each Ti using the adjusted flip angle set for the 
respective T t values. 
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RF transmission field (B/) mapping 

As VFA methods rely on a priori knowledge of the exact flip 
angle in the region of interest, correction for RF transmission field 
(Bj^) non-uniformities is essential. Bj + mapping was done using a 
Bloch-Siegert [32] implementation for phase based Bj + mapping 
[33] using TR—100 ms, Fermi pulse flip angle = 400° (nominal 
Bj =4.35 uT), and a Fermi pulse off-center frequency = 4 kHz. 
Spatial resolution was adapted to the respective Tj mapping 
protocol. The sequence employed double gradient echo acquisi- 
tions with different echo times adapted to the corresponding fat/ 
water frequency shift for each field strength to enable -B 0 -mapping 
[34] . Calculation of transmission field distribution normalized to 
the nominal Bj + was done offline using Bloch simulations in 
Madab that also considered B 0 off-resonances. 

MR hardware 

Measurements were performed on 1.5T/3T/7T whole body 
MRI systems (Avanto/Verio/Magnetom 7T, Siemens Healthcare, 
Erlangen, Germany) running identical software versions. At 1.5T 
and 3T a 12 channel (Siemens Healthcare, Erlangen, Germany) 
head coil was used for signal reception (Rx) while at 7T a 24 
channel Rx coil (Nova Medical, Andover, MA, USA) was 
employed. For transmission (Tx), the integrated body coil was 
used at 1.5T and 3T while the head coil's Tx channel was used at 
7T. 

Validation in Phantom Experiments 

An oval brain like shaped water based phantom (diame- 
ter = 10 cm, length = 20 cm, Tj = 960 ms, T 2 = 100 ms, measured 
at 1.5 T) was built to mimic mean gray and white matter Tj values 
at 1.5T. Manganese chloride was used to adjust T 1; sodium 
chloride was used to adjust the conductivity to that of brain tissue 
[35]. 

To assess the capability of the Bloch simulations to accurately 
emulate slice excitation, the simulated magnetization vector 
profiles were compared to excitation profiles measured exemplar- 
ily at a = 30° at 3.0 T employing a modified FLASH sequence 
with the readout gradient along the slice selection direction. 

Signal intensities over different flip angles calculated from Bloch 
simulations using the actual sequence configuration were com- 
pared to signal intensities derived from standard FLASH 
measurements (TE/TR = 2.5/5 ms; a = 2° — 90°, voxel si- 
ze =(2x2x5) mm 3 ). 

A two-dimensional inversion recovery (IR) prepared technique 
(voxel size: (2x2x5) mm' 5 , non-selective adiabatic hyperbolic 
secant inversion, echo spacing = 5.5 ms, TR = 10 s, a = 90°, turbo 
factor = 5, GRAPPA acceleration factor = 2, 8 inversion times 
of TI = 60/120/240/480/750/1000/3000/5000 ms, scan time = 
16 min) in conjunction with a spin-echo readout served as a 
reference for T! quantification. Ti was calculated offline in Matlab 
using a non-linear least squares three-parameter fit. 

To validate the feasibility and accuracy of the corrected 2D 
VFA T! mapping method accounting for slice profile deforma- 
tions, 2D FLASH images (TE/TR = 2.5/5 ms, scan time = 3.5 s) 
using a set of three adjusted flip angles were acquired. _B/ + -maps 
were used to account for Bj + non-uniformities. 

Ti of the phantom was quantified for all field strengths (1.5T, 
3T, and 7T). For this purpose, a region of interest covering the 
entire cross-section of a central coronal slice of the phantom was 
used. Mean Tj values and standard deviations were calculated for 
this ROI. To elucidate the impact of the signal change due to slice 
profile deformation on the 2D quantification of T^ conventional 
£; + -corrected VFA T! evaluation used for 3D acquisitions that did 
not consider slice profiles were performed using two flip angles 



(2°/ 13°) suggested by the 3D approach [36]. For reasons of 
brevity, the expression "uncorrected 2D VFA" is used for this 
approach in the following. The proposed 2D Bj + corrected VFA 
Ti mapping approach that includes extra slice profile deformation 
versus the conventional approach is called "corrected 2D VFA". 

Ethics Statement 

For the in vivo feasibility study, 6 healthy subjects (mean age: 
28±2 years, 2 females) without any known history of brain disease, 
one MS patient (age: 46 years, female), and one stroke patient (24 
years, male) underwent MR imaging after due approval by the 
ethical committees (registration number DE/CA73/5550/09, 
Landesamt fur Arbeitsschutz, Gesundheitsschutz und technische 
Sicherheit, Berlin, Germany and registration number 
DRKS00003193, 7 Tesla Ultra-High Field Project 7UP, WHO 
International). Informed written consent was obtained from each 
volunteer and patient prior to the study. 

In vivo studies in healthy subjects and patients 

Axial 2D brain images were obtained (FLASH, TE/ 
TR = 2.5 ms/5 ms, scan time = 3.5 s) using flip angle triples, 
which were adjusted to mitigate saturation and noise induced 
errors. Gray matter (GM) and white matter (WM) were segmented 
by thresholding a combination of T! weighted images and T! 
maps derived from the reference Ti measurement. The resulting 
gray matter and white matter masks were applied to the 
uncorrected 2D VFA approach and to the proposed corrected 
2D VFA. The IR-SE based T[ quantification reference method 
described for the phantom measurements was also used for the 
volunteer study. Mean T! and standard deviations were calculated 
for the reference method and for slice profile deformation 
corrected 2D VFA. For volunteer studies, a spatial resolution of 
(2x2x5) mm' 5 was used for data acquisition, which was interpo- 
lated to (1x1x5) mm 3 during reconstruction to improve visual 
appearance. For patient studies, the spatial resolution was adjusted 
to (1.6x1.6x4) mm' 5 during the acquisitions and interpolated to 
(0.8x0.8x4) mm in the reconstruction. 

Results 

Determination of adjusted flip angles using Bloch 
simulations 

Optimum flip angles for T! quantification of the phantom were 
calculated to be 3°/ll o /30° at 1.5T, 3°/ll°/28 0 at 3T, and 3°/ 
ll°/26° at 7T. Adjusted flip angles for T! quantification in vivo 
were found to be 3°/ll 0 /26° at 1.5T, 3°/9°/22° at 3T, and 3°/ 
8°/25° at 7T. 

Monte Carlo simulations 

Figure la shows results derived from Monte Carlo simulation of 
the proposed corrected 2D VFA approach for T] relaxation times 
ranging from 0-3 s. T! precision decreases with increasing Ti, 
while accuracy remains unaffected. If a two flip angle approach is 
used as a reference, the T! precision index was improved for the 
three flip angles approach by 9% for T! values ranging from 
100 ms to 300 ms, by 13% for T! values ranging from 300 ms to 
1000 ms, and by 12% for T! values ranging from 1000 ms to 
3000 ms. 

Validation in phantom experiments 

Bloch simulated and measured excitation slice profiles derived 
from the phantom are depicted in Figure lb for an excitation pulse 
of a = 30°. Simulation and measurement agree very well. The 



PLOS ONE | www.plosone.org 



3 



March 2014 | Volume 9 | Issue 3 | e91318 



Rapid 2D T1 Mapping at 1.5T, 3T, and 7T 




Figure 1. Monte Carlo simulation and impact of slice profile deformation on FLASH signal, (a) Simulated density distribution 
(logarithmic) of true T, values versus estimated T q values derived from corrected 2D VFA using three flip angles. For each T, value, 1 0000 experiments 
were performed. The use of three flip angles improved the T, precision index by 1 1 % compared using only two flip angles, (b) Bloch simulated (blue 
solid line) and measured (black dashed line) slice profiles (signal relevant component) in the steady state of a FLASH sequence show good agreement. 
The slice profile substantially deviates from an ideal rectangular slice profile, (c) Theoretical (green dashed line), Bloch simulated (blue solid line), and 
measured (black dots) signal intensities for FLASH as a function of the flip angle (at T q = 1 000 ms). The Bloch simulated signal intensities agree very 
well with the phantom measurements. 
doi:1 0.1 371 /journal.pone.0091 31 8.g001 



saturation dependent slice profile deformation is characterized by 
the RF pulse shape, TR, T] and the flip angle. The slice profile 
deformation and B t + non-uniformity induced discrepancy between 
theoretical and measured signal curves for FLASH is demonstrat- 
ed in Figure lc. While theoretical signal intensities using Eq. 2 
deviated substantially from the measurements, the signal curve 
derived from the B] + corrected simulations including slice profile 
considerations perfectly matched the experimental data. The flip 
angle providing maximum signal was found to be 11° (considering 
slice profile deformations), which is almost twice the theoretical 
Ernst angle of 6°. 

Bj maps of the phantom at 1 .5T, 3T, and 7T are shown in 
Figure 2a. As expected, the Bj + homogeneity decreases with 
increasing field strength. At 1.5T, T[ quantification using the 
inversion recovery approach as a reference yielded 
T 1 = (963±7)ms for the cross-sectional ROI of the phantom. At 
3T, the reference Tj value was found to be ( 1 022 ± 10)ms. For the 
same phantom, T q = (1002± 10)ms was obtained at 7T as 
illustrated in Figure 2b. Conventional, only Bj + corrected 2D 
VFA T q mapping using Eq. 2 substantially underestimated T| and 
yielded Ti = (372±28)ms at 1.5T, Tj = (376±23)ms at 3T, and 
T q = (357±60)ms at 7T as demonstrated in Figure 2b. The 
inhomogeneous T q map at 7T also underlines that a solitary Bj + 
correction is not sufficient to compensate for the non-linear flip 
angle dependent signal changes. In comparison, slice profile 
deformation and Bj + corrected 2D VFA provided 
T 1 = (966±67)ms at 1.5T, T q = (1010±47)ms at 3T and 
T q = (1017±89)ms at 7T as illustrated in Figure 2b, which 
matches the reference data fairly well. The proposed slice profile 
corrected 2D VFA approach yielded a mean deviation of T] from 
the reference value as low as (0.2 ±1.4)% averaged over all field 
strengths (Figure 2c). In comparison, an error in T[ quantification 
of (— 63 ± 1.5)% was observed for the 2D VFA technique without 
slice profile deformation correction as shown in Figure 2c. 

All measurements, including Tj reference measurements, 
showed inconsistent results in the neck of the water bottle at 7T, 
because this region was outside of the covered volume of the local 
transmission coil. 



Volunteer and Patient Studies 

T] mapping performed in healthy volunteers at 1.5T, 3T, and 
7T using the IR-SE reference measurement, the uncorrected 2D 
VFA, and the corrected 2D VFA are depicted in Figure 3. T q 
values for gray matter and white matter are illustrated in Table 1 . 

At 1.5T, flip angles were 3°/l l°/26° for the proposed corrected 
2D VFA and 2°/ 13° for uncorrected 2D VFA. The absolute mean 
difference between IR reference and uncorrected 2D VFA was 
60% for gray matter and 56% for white matter. The absolute 
mean difference between IR reference and corrected 2D VFA was 
1 1 % for gray matter and 1 7 % for white matter. 

At 3T, flip angles were 2°/ 12° for uncorrected 2D VFA and 3°/ 
9°/22° for the proposed corrected 2D VFA. The absolute mean 
difference between IR reference and uncorrected 2D VFA was 
68% for gray matter and 59% for white matter. The absolute 
mean difference between IR reference and corrected 2D VFA was 
1 1% for gray matter and 6% for white matter. 

At 7T, flip angles were 2°/l 1° for uncorrected 2D VFA and 3°/ 
8°/25° for the proposed corrected 2D VFA. The absolute mean 
difference between IR reference and uncorrected 2D VFA was 
61% for gray matter and 57% for white matter. The absolute 
mean difference between IR reference and corrected 2D VFA was 
22% for gray matter and 44% for white matter. 

The acquisition time for the IR-SE reference measurement was 
16 minutes for each slice at each field strength. For corrected 2D 
VFA T[ mapping, acquisition time including Bj + mapping was 
18 s for each slice. 

The Tj measurement of the MS patient at 7T revealed 
(2758±615)ms in three periventricular lesions vs. (1842±62)ms in 
the surrounding white matter (Fig. 4a). T q in the subcortical lesion 
of the stroke patient was (2554±316)ms vs. (951±74)ms of the 
surrounding white matter at 3T (Fig. 4b). 

Discussion 

This study demonstrates the feasibility of rapid 2D T t 
quantification using variable flip angles in phantom experiments 
and in vivo at magnetic field strengths of 1.5T, 3T and 7T. For this 
purpose, implications of RF transmission field non-uniformities, as 
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Figure 2. Normalized B 7 * maps and T, maps of the water 
phantom at 1.5T, 3T, and 7T. (a) Measured B* normalized to the 
nominal 6, + (4,35 uT) in percent, (b, top row) T q derived from the 
inversion recovery reference method (T q = (963±7)ms at 1.5T, 
T, =(1022±10)ms at 3T, T, = (1002±1 0)ms at 7T). (b, center row) 
Uncorrected 2D VFA approach that does only consider B, + but not slice 
profile deformations. Adjusted flip angles (a, =2°, ot 2 = 13°) were used. 
This approach revealed a T, underestimation of -63.0±1.5% averaged 
over all field strengths (T, = (372±28)ms at 1 .5T, T, = (376±23)ms at 3T, 
T q = (357±60)ms at 7T) compared to the reference measurement, (b, 
bottom row) With the proposed 6, + corrected and slice profile 
distortion corrected 2D VFA measured a T, deviation of only 0.2±1.4% 
averaged over all field strengths (T q = (966±67)ms at 1.5T, 
T, =(1010±47)ms at 3T, T, =(1017±89)ms at 7T) compared to the 
reference measurement, (c, top row) Difference map between 
uncorrected 2D VFA and IR reference, (c, bottom row) Difference 
map between corrected 2D VFA and IR reference. 
doi:1 0.1 371 /journal.pone.0091 31 8.g002 

well as excitation slice profile deformations due to non-ideal RF 
pulses were carefully considered. For correction of slice profile 
deformation induced errors, Bloch simulations were employed. 
Corrected 2D VFA T, mapping was found to be accurate when 
compared to reference inversion recovery acquisitions in phantom 
studies at all tested field strengths. For T t mapping of the brain, 
mean T[ values measured by corrected 2D VFA show a fairly 
good agreement with the reference measurement at 1.5T and 3T 
for both gray and white matter. T] relaxation times found for 
corrected 2D VFA at 7T were larger than that derived from IR- 
SE. The applicability and efficacy of the proposed procedure to 
calculate adjusted flip angles in order to mitigate noise induced 
errors was confirmed by our simulations and experimental data. 

The inclusion of Bj + maps made 2D VFA T q mapping immune 
to non-uniformities in the Bj + transmission field, which are 



normalized B 




200 

100 <s ~ 
50 I" 



B 

IR reference 

uncorrected 
2D VFA 




• 


• 


H 






corrected 
2D VFA 


'■ . ' 








$ 





J 



1.5 



0.5 



difference 
IR reference- 
uncorrected 
2D VFA 



difference 
IR reference- 
corrected 
2D VFA 




0 ~ 



Figure 3. B, + mapping and T, mapping in healthy volunteers at 
1.5T, 3T, and 7T. Exemplary normalized 6, + maps (a) and T q maps (b) 
of a healthy brain using the inversion recovery reference measurement, 
the uncorrected 2D VFA, and the proposed corrected 2D VFA approach 
at 1.5T, 3T, and 7T. The difference maps between IR reference and 
uncorrected/corrected 2D VFA are illustrated in (c). 
doi:10.1371/journal.pone.0091318.g003 

pronounced at high and ultrahigh field strengths. Magnitude 
based Bj + mapping methods such as double angle [37,38,39,40] or 
the use of inversion recovery prepared extra scans [21,41] have 
been presented previously to correct for flip angle inaccuracies. 
However, most of these techniques rely on equilibrium magneti- 
zation before each acquisition or need extra calibration routines 
and therefore add extra scan time. Flash-EPI hybrid sequences 
employing adiabatic pulses [42] were proposed to potentially 
provide means for overcoming B 1 inhomogeneities while 
maintaining reasonable scan time. However, SAR constraints 
dictate flip angle limits and adiabatic pulses demand longer echo 
times especially at higher magnetic fields, which lead to partial 
dephasing of the magnetization. The phase based Bj + mapping 
technique used here is fast and reliable, even in low signal regions 
[33] of the brain. The rapid data acquisition reduces motion 
sensitivity in Bj mapping, which is ultimately beneficial for in vivo 
T! quantification. 

The reference method was already accelerated by using parallel 
imaging with an effective reduction factor of R = 2 together with a 
spin echo turbo factor of 5 and a reduced phase field of view. The 
corrected 2D VFA implementation used for the volunteers offers a 
scan time advantage factor of approximately 570 versus the 
unaccelerated inversion recovery reference method. Obviously, 
corrected 2D VFA supports modest accelerations, which would 
further enhance the speed advantage over the reference. 
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Table 1. T, values for gray matter and white matter in healthy volunteers at 1.5T, 3T, and 7T derived from IR-SE reference 
measurement, uncorrected 2D VFA, and corrected 2D VFA. 





T, [ms] 






1.5T 


3T 


7T 


white matter 


IR reference 


678±10 


91 1 ± 1 5 


1284±22 




uncorrected 2D VFA 


298±14 


377±25 


554±21 




corrected 2D VFA 


791 ±21 


969 ±85 


1855+141 


gray matter 


IR reference 


1154±82 


1615±149 


2065 ±69 




uncorrected 2D VFA 


465 ±40 


518±42 


804 ±48 




corrected 2D VFA 


1282 ±78 


1433±80 


2524±137 



doi:1 0.1 371/journal.pone.0091 31 8.t001 



Notwithstanding the utility of the reference method that relied on 
equilibrium magnetization before each inversion pulse for T! 
quantification, the potential to shorten scan times is limited. As T! 
of brain tissue increases with field strength, the equilibrium 
condition of TR>5xTj requires even longer TRs, which further 
increase scan time. Unlike the reference method, moving to higher 
magnetic fields does not add a scan time penalty to corrected 2D 
VFA. Even more, the lack of inversion pulses in 2D VFA reduces 
magnetization transfer effects and RF power deposition. 

Our results demonstrate severe T] quantification errors due to 
slice profile deformations in the uncorrected 2D VFA Tj mapping 
using commonly accepted signal equations. This underscores that 
theoretical signal predictions only hold true for ideal rectangular 
slice profiles assuming uniform excitation over the entire slice or 
volume. Slice selection RF pulses are however bandwidth limited 
and lead to a spatially non-uniform flip angle distribution in slice 
selection direction. This evokes T r , TR-, and flip angle dependent 
saturation phenomena that create distorted slice profiles. Conse- 
quently, it is essential to correct for slice deformations before T! 
values derived from 2D VFA acquisitions can be considered 
accurate. 

The match of T] values derived from IR-SE reference 
measurements with T! values derived from the corrected 2D 
VFA at all tested field strengths suggests that the sequence related 
deviations were successfully addressed by our 2D VFA implemen- 
tation. In-vivo T! maps obtained with the two methods, however, 
revealed some residual differences. This underscores the challeng- 
es of parametric mapping when moving from phantoms to in vivo, 



where tissue characteristics become relevant and play a role for T! 
quantification. It is to be expected that magnetization transfer 
weighting, weighting of T! relaxation components, and T 2 
relaxation components in mixed tissue is different for different 
Ti mapping techniques, which ultimately leads to different mean 
T! values. Therefore, it is challenging to establish normal values 
for tissue that are valid across all Tj quantification techniques. The 
broad variety of reference Tt values found in the literature 
supports this statement [25,26,27,28,29,30,31]. We also performed 
simulations using a three pool model of mixed white matter tissue 
(myelin, myelinated axons, and mixed water pools) together with 
volume fractions reported for adults [43] . The difference in mean 
T t between the IR-SE reference and the proposed corrected 2D 
VFA was found to be 15%. In our measurements, the linear trend 
of increasing T [ with increasing field strength could not be shown 
for the corrected 2D VFA method at 7T. However, for all field 
strength the discrepancy between the reference measurement and 
corrected 2D VFA T! mapping (mean deviation of 18%) was 
reduced versus the discrepancy between the reference and 
uncorrected 2D VFA T[ mapping (mean deviation of 63%). 

Admittedly, it may take several excitations before the magne- 
tization reaches steady state in gradient echo techniques, which 
depends on the RF pulse shape, T! relaxation time, TR, and the 
flip angle used. For this reason, we integrated 200 dummy pulses 
for T! mapping to approximate steady state. This approach 
resulted in total scan times of approximately 1 8 s for single slice 
2D VFA T! mapping in the phantom and in healthy volunteers 
(spatial resolution = 2 x2 x5 mm ) including all three flip angles 
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Figure 4. Patient measurements using the corrected 2D VFA method, (a) T q -map of a multiple sclerosis patient at 7T using the corrected 2D 
VFA method. The map shows three periventricular lesions (white arrows). Mean T q of the lesions was (2758±615)ms vs. (1842±62)ms in the 
surrounding white matter, (b) Corrected 2D VFA T,-map of a stroke patient derived from 3T acquisitions showing a subcortical lesion (white arrow). 
The stroke occurred five months before measurement. T q of the lesion was (2554±316)ms vs. (951 ±74)ms in the surrounding white matter. 
doi:10.1371/journal.pone.0091318.g004 
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and including mapping. For the patient measurements, a 
clinically acceptable scan time of 30 s was achieved for single slice 
T] mapping (spatial resolution = 1.6x1.6x4 mm' 5 ). As all variable 
flip angle Ti quantification techniques rely on the steady state 
approximation, fast moving blood might not allow the steady state 
to form, alter the signal, and thus bear the potential to deem blood 
Tj quantification inaccurate. 

In order to maintain consistent scan protocols throughout all 
field strengths, we used an echo time of 2.5 ms and a repetition 
time of 5 ms due to specific absorption rate (SAR) considerations 
at 7T. It is a recognized limitation that this approach is suboptimal 
at 1 .5T and 3T. Here, shortening of TR would further reduce 
scan time. Moreover, RF Pulse shapes, gradient strengths and slew 
rates, receiver bandwidth, and spatial resolution can be adapted to 
use the full range of the system's specification at 1.5T, 3T and 7T. 
The repetition times and scan times reported here for B 0 and B 1 
mapping were limited by SAR constraints at 7T and of course can 
be further shortened at 1.5T and 3T. Implementation of these 
changes affords mapping for a single slice in 5 s including B t + 
mapping. To reduce the computational effort, an accelerated 
algorithm for T, quantification was implemented in Madab, 
which enables generation of T[ maps using a matrix size of 
256x256 pixels in less than three seconds. 

The proposed approach for the calculation of adjusted flip 
angles to reduce T! quantification errors in 2D VFA T! mapping 
by balancing dynamic range and SNR does not yet account for 
Bj + inhomogeneities or RF transmitter mis-calibrations, which can 
be hardware-, setup-, subject-, and also organ specific. Therefore, 
the proposed approach requires proper RF transmitter adjustment 
- a procedure which is included in the pre-scan calibration of 
clinical scanners - in order to get as close as possible to the adjusted 
flip angles. The choice of flip angles depends on the grade of slice 
profile deformation and therefore varies with the RF pulse shape. 
Our simulations and experiments focused on a bandwidth 
truncated sine RF pulse with a bandwidth of 1 kHz and a 
duration of 2 ms. For other RF excitation pulses, adjusted flip 
angles might be different and require a onetime pre-calculation 
before used for T] mapping. In contrast to the theoretical 
assumptions based on signal equations, the deformation of the slice 
profile in 2D gradient echo acquisitions shifts the suggested flip 
angles to larger values. Consequently, these flip angles may be 
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